UNCLASSIFIED 

u 429178 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


'4ll||l||i|HIIINI|l|lllllll|lllNINHlHIH|IIIUIilMIMNIIIIilHIIIW|l|HIHINIIINIII|l|IM|<i  -llrilIlillHIIIHIIIIIIIiIIIIHfiilnmmHli;.illllNluMyiHlllllllHlilMimumtiiMHiiiiiiim«ii;iii.iiHiiiiiii . . . 


J0TIC1:  Whan  fovenawnt  or  other  drevinia,  eneci- 
ficatlone  or  other  data  are  uaed  for  toy  purpoee 
other  than  In  connection  vlth  a  definitely  related 
goveraaent  procuiwent  operation,  the  U.  8. 
Qoyejment  thereby  Incur*  no  mpooaibility,  nor  any 
oblimtion  ’nhat.oever;  and  ttaa  fact  that  thi  Oovern- 
aent  nay  have  fonulated,  fumlahed,  or  In  any  vav 

J””1.1"1 ipwltlartleM,  or  otter 
li  not  to  te  rofudod  tar  lnpllcotloo  or  otter- 
elae  ae  in  any  aamaer  licensing  the  bolder  or  any 
other  peraon  or  corporation,  or  conveylnc  any  riahta 
or  peraiasion  to  ninufacture,  use  or  aell  any 
patented  invention  that  nay  in  any  eey  he  related 
thereto. 


O* 

CJ> 

Q 


03 

o 

•  u 

CO 

o 


CO 


01 

£> 


IS.**-/ 


TECHNICAL  INFORMATION  SERIES 


R63SD58 


SOME  ASPECTS  OF  TURBULENT  SCATTERING 
OF  ELECTROMAGNETIC  WAVES  BY 
HYPERSONIC  WAKE  FLOWS 


DDC 

J  F.fcb  1 7  Wb4 

buoi-’J  -2  *■ 

TU»iA  D 


fj 

i 


K.T.  YEN 


SPACE  SCIENCES  LABORATORY 

GENERAL  @  ELECTRIC 

MISSILE  AND  SPACE  DIVISION 


CONTENTS 


I  . 

ABSTRACT 

I.  INTRODUCTION 

II.  SCATTERING  OF  ELECTROMAGNETIC  WAVES  BY 
TURBULENT  MEDIUM 

III.  SOME  ASPECTS  OF  THE  WAKE  TURBULENCE 
PHENOMENON 

IV.  SCATTERING  CROSS  SECTION 

V.  CONCLUSIONS  AND  RECOMMENDATION  FOR 
FURTHER  WORK 

APPENDIX  I 

APPENDIX  II 

REFERENCES 

FIGURES 


ACKNOWLEDGEMENT 


ABSTRACT 


The  present  analytical  study  of  the  turbulent  scattering  of  electro¬ 
magnetic  waves  is  directed  to  certain  features  of  the  phenomenon  considered 
as  peculiar  to  hypersonic  wakes.  Effects  of  flow  intermittency,  non-isotropy 
of  the  turbulence  structure,  and  finite  width  of  the  wake  on  the  scattering 
cross  section  in  its  frequency  and  aspect-angle  dependence  are  of  primary 
concern. 

Phenomena  of  turbulence  revelant  to  the  scattering  problem  are 
briefly  discussed.  In  particular,  it  is  shown  that  electron-density  fluctua¬ 
tions  of  a  turbulent  nature,  in  addition  to  those  caused  by  "turbulence", 
will  be  produced  by  the  intermittency  behavior  of  turbulent  wake  flows. 

This  fluctuation  is  found  to  depend  on  the  mean  electron-density  distribu¬ 
tion,  and  Townsend's  intermittency  "ft"  function. 

Phenomenological  consideration  of  the  physical  process  by  which 
the  turbulent  electron-density  fluctuations  are  likely  to  be  produced  in 
hypersonic  wakes  yields  an  expression  for  the  intensity  of  the  fluctuations 
in  terms  of  the  gradient  of  the  mean  electron -density  distribution,  the 
intensity  of  the  turbulent  velocity  fluctuation,  and  a  turbulence  scale  of  the 
turbulent  velocity  field. 

Application  of  the  above  considerations  to  turbulent  scattering  by 
undcrdense  wakes  has  been  made.  In  contrast  to  the  ionospheric 
scattering,  two  contributions  to  the  scattering  cross  section  are  obtained 
in  the  present  case:  the  first  one  arises  from  the  intermittency 
phenomenon  and  vanishes  if  the  intermittency  is  omitted,  and  the  second 
one  is  due  to  "turbulent"  fluctuations.  The  frequency-dependence  of 
these  two  contributions  are  seen  to  be  quite  different.  It  is  further  shown 
that  flow  intermittency,  non-isotropy  of  turbulence  structure,  and  finite 
width  of  the  wake  will  introduce  aspect-angle  dependence  into  the  scatter¬ 
ing  cross  section.  A  possible  form  of  the  total  scattering  cross  section 
for  hypersonic  wakes  is  also  given. 

Much  of  the  problem  remains  to  be  studied.  Areas  of  additional 
work  of  immediate  concern  are  briefly  described. 
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I.  INTRODUCTION 

Interest  in  turbulence  as  a  mechanism  for  scattering  of  electro- 
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magnetic  waves  probably  began  when  Booker  and  Gorden  suggested  it 
to  explain  some  anomalous  phenomena  of  radio  scattering  in  the 
troposphere.  The  original  theory  developed  by  Booker  and  Gorden 

A  ij  j 

formed  the  basis  for  further  theoretical  and  experimental  work  *  ’  ’  in 

many  related  scattering  problems.  As  applied  to  ionosphere  back- 
.  3 

scattering  ,  the  essence  of  this  theory  is:  when  the  scattering  volume 
is  much  larger  than  Xq  ,  where  Xq  is  the  wave  length  of  the  incoming 
electromagnetic  wave,  and  when  the  electron  density  is  such  that  the 
plasma  frequency  is  less  than  the  wave  frequency,  (i.  e.  the  plasma  is 
underdense),  turbulent  irregularities  in  electron  density  responsible 
for  the  back-scattering  are  those  with  wave  vector  equal  to  2k,  where 
k  is  the  wave  vector  of  the  incident  electromagnetic  wave. 

In  determining  the  scattering  cross-section,  however,  assumptions 
concerning  the  intensity  of  the  electron  density  fluctuations  and  the  form 
of  the  correlation  or  spectrum  function  must  be  made  on  account  of 
insufficient  knowledge  about  the  turbulence.  Several  different  formulas 
have  been  proposed  for  the  scattering  cross  section  by  different  workers, 
but  none  of  these  appear  to  be  entirely  satisfactory.’’  Hence,  at  the 
present  time  the  theory  of  turbulent  scattering  of  radio  waves  in  the 
ionosphere  must  be  considered  as  still  incomplete. 

It  is  also  recognized  that  turbulent  scattering  may  be  responsible 
for  radar  return  signals  from  the  wakes  of  high-speed  reentry  objects 
at  certain  altitudes.  Experimentally,  the  back-scattering  signals  from 
the  wake  of  a  reentry  object  have  been  found  to  be  undetectable  by  radar 
until  the  object  is  at  an  altitude  at  which  the  wake  apparently  becomes 
turbulent.  It:  is  reasonable  to  expect  no  significant  back-scattering 
returns  from  a  laminar  wake  at  small  aspect  angles  with  respect  to  the 
pi  opagation  direction  of  the  radar  wave  (Figure  1),  (this  problem  was 
investigated  in  reference  6  based  on  a  simplified  model). 


1 


. . Ill . I . 'I'""1"" . "I"' . . . 


The  problem  of  turbulent  scattering  by  underdense  wakes  has  been 

studied  in  references  7,  8,  9  and  others.  The  theoretical  basis  of  these 
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studies  can  be  traced  back  to  the  work  of  Booker  and  Gorden  ,  and  does 

3 

not  appear  to  have  progressed  much  beyond  Booker's  recent  work  on 
ionspheric  radio  propagation.  The  purpose  of  this  study  is  to  examine 
the  physical  nature  of  the  scattering  problem  peculiar  to  the  wakes  of 
reentry  objects  of  certain  type.  A  new,  but  tentative,  form  of  the 
scattering  cross  section  has  been  obtained.  Since  solution  to  this 
problem  is  not  yet  complete,  this  report  may  serve  to  define  the  scope 
and  present  status  of  this  problem,  and  also  areas  of  future  work. 

II.  SCATTERING  OF  ELECTROMAGNETIC  WAVES  BY  TURBULENT 
MEDIUM 

In  the  theory  originated  by  Booker  and  Gordon*,  irregular  fluctua¬ 
tions  of  the  dielectric  constant  of  a  medium  are  considred  as  the  sources 
for  scattering  of  electromagnetic  waves  by  the  medium.  As  in  the 
turbulence  theory***,  the  dielectric  constant  of  a  medium  may  be 
separated  into  its  mean  value  K  and  its  flucation  AK,  where  both  K  and 
AK  may  be  functions  of  the  space  coordinates  r  (or  x.,  i  =  1,  2,  3)  and 
the  time  t. 

When  a  plane  trasverse  electromagnetic  wave  of  the  form 
A()  exp  {  -  i(wt  -  k  •  r )  )  is  incident  on  a  volume  V  of  such  turbulent 
medium,  the  scattering  cross  section  (per  unit  solid  angle)  can  be 
written  in  the  following  form 


4  2 

k  sin  $ 

16  tr2 


AK  >  AK 

v-r(ri)ir 


V  V 


i(k  -  km)  •  (f  -  ?  ) 

(rje  dV.  d V 

1  z 


(1) 


where  m  =  (r/r)  is  a  unit  vector  directed  from  the  origin  of  coordinates 
(chosen  within  the  scattering  volume  V )  to  the  observation  point  (Figure 
2),  |8  is  the  angle  between  ,Aq  and  m,  and  the  mean  value  of  a  quantity 
is  indicated  by  the  symbol  <  >  which  yields,  in  the  present  case,  the 
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correlation  function  for  the  dielectric -constant  fluctuations.  The  above 

4 

formula  can  be  easily  deduced  from  Tatar  ski's  results. 

In  obtaining  Eq.  (1),  a  perturbation  technique,  in  principle  the  same 
as  "Born  approximation"  in  its  first  approximation,  has  been  used  to 
solve  the  Maxwell  equations.  As  a  consequence,  Eq.  (1)  becomes 
invalid  if  K  vanishes  in  the  scattering  volume.  Let  the  dielectric  con¬ 
stant  be  of  the  form 


K 

TT  = 

v 


(2) 


where  K  denotes  the  value  ofK  in  vacuum,  and  CO  is  the  plasma 
v  p 

frequency.  Thus,  Eq.  (1)  applies  only  to  "underdense"  plasmas,  i.  e. 

(0  <<  (0  .  In  addition,  since  both  A  and  k  are  taken  as  constant  in  the 
P  o 

scattering  volume,  the  quantity  grad  f  in  K)  must  remain  small.  However, 
the  fluctuation  in  electron  density  An  is  allowed  to  be  comparable  to  N, 
the  mean  electron  density.  This  is  due  to  the  fact  that 


AK 

K 

v 


(3) 


a  2  2 

and  will  remain  small  compared  to  1  for  An  —  N  as  long  as  (0  <<  to  . 

P 

(Large  mass -density  fluctuations  have  been  observed  experimentally,  see 

2 

reference  11.  Speculations  have  been  put  forward  that  <  (An)  >  may  be 
comparable  to  N  ) 

For  back-scattering,  i.e.  m  =  -  k,  and  j9  =  (l/2)ff,  the  scattering 
cross  section  reduces  to 


0  = 


I6ir 


.  ak  ak  ^ 

K  K  r2  ' 


i  2k 


r'i  -v 


V  V 


dVl  dV2 


(4) 


As  indicated  by  Tatarski  ,  a  change  of  variables  can  be  introduced  as 
follows 
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Then 
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where 


R  = 


+  r 


P  *  r  j 


i  2  k  •  P 

B(R  ,  P )  e  dV 

P 


B{R,  P)  »  <  Trj)  > 
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Under  the  conditions  that  ft)  <<  u  and  the  fluctuation  in  dielectric 

r 

constant  is  due  to  that  of  the  electron  density,  as  given  by  Eq.  (  3), 
Eq.  (5)  becomes 


(7  3 


i  2  k  •  *p 

An(R)  An(R  +  P)  >  e  dVp 


(6) 


where  r  {  =  e  /c  m  )  is  the  classical  electron  radius. 

If  the  scattering  volume  V  is  very  large,  so  that  the  integral  in 
Eq.  (5)  with  respect  to  P  can  be  considered  as  Fourier  integral,  the 
following  result  is  obtained:  The  contribution  to  the  back -scatter  in  a 
comes  from  one  (and  only  one)  Fourier  component  of  the  fluctuations  of 
the  dielectric  constant,  i.e.  that  component  whose  wave  vector  equals 
to  twice  the  wave  vector  of  the  incident  electromagnetic  wave.  This 
result  is  generally  considered  as  applicable  to  scattering  problems  for 
which 


L  > 
n 


(7) 


where  X  is  the  length  of  the  incoming  electromagnetic  wave,  L  the 

°  '  n 

scale  of  the  largest  turbulent  "eddies”  of  the  dielectric  constant  fluctua- 


4 


tions,  and  I  that:  of  the  smallest  turbulent  "eddies"  (these  "eddies"  are 
n 

not  meant  for  flow  eddies).  The  size  of  the  scattering  Volume  V  must 

3 

be  necessarily  larger  than  Lr  . 

Determination  of  the  correlation  function  is  known  to  be  an  out¬ 
standing  problem  in  the  turbulence  theory.  A  discussion  of  this  problem 
will  be  presented  in  the  next  section.  In  the  case  when  the  turbulence  is 
homogeneous,  i.e.  when  the  correlation  function  B  depends  only  on  (0, 
the  expression  for  (7,  Eq.  (6),  reduces  to 


i  2  k  •  P 


B(0)  e 


dV 


where  B  (P)  *  <  An(R)  AnfR  +  P)  >  depends  only  on  P  . 

If  the  fluctuations  An(R)  and  An  (R  +  p)  are  not  correlated,  then 

<  An(R)  An(R  +  P)  >  will  be  zero,  except  when  P  ■  0.  In  this  case, 

2 

the  back-scattering  cross  section  (for  a  perfect  gas,  since  <(An)  t>  is 
equal  to  N/V  )  becomes 


a  =  r  2  N  V  .  (8) 

l  e  ’  ' 

12 

This  is  the  cross  section  for  "incoherent"  back-scattering.  It  is 
generally  recognized  that  turbulent  scattering  is  far  from  being 
incoherent,  since  the  fluctuations  are  not  quite  random.  It  seems 

reasonable  to  suggest  the  turbulent  scattering  cross  section  to  be  the 

2  2  2  2 
sum  of  an  incoherent:  part  ar  N  V  and  a  coherent  part  br  N  V 

where  a  and  b  are  presumbly  independent  of  the  mean  electron  density 

N.  It:  will  be  of  interest  to  see  whether  this  form  of  turbulent  cross 

section  can  be  justified  with  our  present  incomplete  understanding  of 

h o m o g e n e o u s  turbulence. 

HI.  SOME  ASPECTS  OF  THE  WAKE  TURBULENCE  PHENOMENON 
A  straightforward  way  of  solving  the  wake  scattering  problem  is 
to  assume  a  form  for  the  correlation  function  in  evaluating  the  integral 
in  Eq.  (5).  While  the  correlation  coefficient  can  usually  be  approxi- 
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mated  by  suitable  formulas  valid  under  certain  circumstances,  to 
predict  the  intensity  of  the  electron -density  fluctuations  presents  a  more 
difficult  problem.  In  the  wake  of  a  hypersonic  reentry  body,  the  physical 
processes  by  means  of  which  the  electron-density  fluctuations  are  pro¬ 
duced  appear  to  be  sufficiently  different  from  those  present  in  the  iono¬ 
sphere  so  that  an  examination  of  them  is  warranted. 

First,  however,  the  distinction  between  a  turbulent  velocity  field 
and  a  turbulent  scalar  field  should  be  noted  (see  reference  10  for  the 
isotropic  case).  The  correlation  function  for  a  velocity  field  is  a 
second -order  tensor,  while  that  for  a  scalar  field  is  a  scalar.  The 
physical  processes  of  "dissipation"  for  these  two  Helds  also  differ 
significantly.  Therefore,  it  is  important  to  recognize  the  difference 
in  form  between  their  correlation  or  spectrum  functions  (this  point  was 
over  looked  by  several  workers).  The  turbulence  scales  of  the  scalar 
field  may  also  differ  from  those  of  the  velocity  field. 

As  generally  recognized,  significant  back-scattering  by  hypersonic 
wakes  apparently  takes  place  only  when  the  wake  flow  is  turbulent.  Thus, 
the  fluctuations  of  electron  density  are  caused  by  the  turbulent  motions 
of  the  velocity  field,  i.  e.  ,  the  action  of  the  turbulent  eddy  motions  on 
the  non-uniform  distribution  of  electron  density  is  considered  as  the 
mechanism  for  "turbulent  fluctuations"  of  electron  density.  By  ele¬ 
mentary  phenomenological  consideration,  an  expression  for  the  fluctua¬ 
tion  in  terms  of  the  local  electron -density  gradient,  and  the  scale  and 
intensity  of  the  turbulent  motion  will  be  given.  It  may  be  mentioned  that, 
in  the  theory  of  incoherent  scattering,  formulas  for  the  electron-density 
fluctuations  have  been  worked  out  (see  references  12,  13  and  14)  by 
means  of  more  rigorous  methods. 

In  the  wake  of  a  high-speed  reentry  body,  the  turbulent  flow  is 
known  to  be  non-isotropic  and  non-homogeneous.  The  turbulent  scalar 
field  is  also  expected  to  be  non-isotropic  and  non-homogeneous,  i.  e. 
the  turbulent  scales  in  the  longitudinal  and  transverse  directions  are 
different,  and  the  turbulent  intensity  varies  over  the  wake  (decaying  on 
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15 

account  of  dissipation).  It  is  generally  accepted  ,  however,  that  when 
the  turbulent  Reynolds  numer  is  high  enough,  there  exists  a  range  of 
small  eddies  (high  wave -numbers)  for  which  the  turbulence  is  essentially 
isotropic.  If  the  eddies  responsible  for  the  scattering  are  in  this  range, 
the  theory  of  locally  isotropic  turbulence  may  be  used. 

When  the  length  of  the  incoming  electromagnetic  wave  is  com¬ 
parable  to  smaller  than  the  wake  diameter,  however,  the  larger  eddies 
become  the  chief  sources  of  scattering.  The  effect  of  non-isotropy  may 
become  significant,  e.  g. ,  the  scattering  cross  section  may  become 
aspect-angle  dependent.  Moreover,  the  condition  given  by  Eq,  (7)  is 
no  longer  satisfied,  and  the  volume  of  integration  Eq.  (6)  will  be  a 
finite  one.  Consequently,  finite -width  effect  of  the  wake  should  appear 
in  the  scattering  cross  section,  with  the  ratio  between  the  wave  length 
and  the  wake  diameter  as  a  parameter.  ThiB  effect  will  be  treated  in 
the  next  section. 

A  prominent  feature  of  turbulent  wake  flows  is  the  phenomenon  of 
1 6 

intermittency,  It  will  be  shown  in  the  following  that  electron -density 
fluctuations,  in  addition  to  those  caused  by  the  "turbulence",  will  be 
introduced  by  the  intermittency. 

(1)  Effect:  of  Intermittency  on  Electron  Density  Fluctuation 

According  to  Townsend  ,  the  cause  of  the  flow  intermittency  is 
the  production  of  a  convoluted  boundary  surface  between  the  turbulent 
and  non -turbulent  fluid  by  the  large  eddies  (of  fluid  motion).  Thus, 
across  the  wake  an  intermittency  factor  y  is  introduced  to  show 


statistically  the  fraction  of  the  time  that  the  flow  is  found  to  be  turbulent. 

Since  the  hot  boundary -layer  gas  goes  into  the  turbulent  core  of 

the  wake,  electrons  produced  by  thermal  ionization  will  be  presented  in 

the  turbulent  core.  In  the  non-turbulent  fluid,  the  electron  density  is 

* 

very  low  and  may  be  taken  as  zero.  Intermittency  of  the  turbulent:  flow 


>1* 

Although  this  condition  appears  to  be  valid  for  slender  bodies,  in  general 
enthalpy  in  the  outer  inviscid  wake  is  not  negligible  comparable  to  the 
inner  turbulent  wake  for  hypersonic  blunt  bodies17.  The  derivation  and 
results  to  follow  can  be  readily  modified  to  include  the  electron  density 
distribution  in  the  non-turublent  fluid. 
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will,  consequently,  introduce  at  a  fixed  point  fluctuations  of  electron 

density,  in  addition  to  those  present  in  the  entirely  turbulent  fluid, 

Let  6  (r ,  t)  denote  the  probability  of  turbulence  at  r  ,  such  that 

ft  is  aero  in  non-turbulent  fluid  and  is  unity  in  turbulent  fluid,  ^  Time 

average  of  6  yields  y(r ),  the  intermittency  factor  at  r .  Let  n(r ,  t) 

be  the  electron  density  in  the  turbulent  fluid,  N(r  )  the  mean  value  of  n 

averaged  {in  time)  over  both  the  turbulent  and  the  non* turbulent  fluid, 
w  1  / 
and  n(r  )  the  mean  value  of  n  in  the  turbulent  fluid,  Townsend  has 

shown  that  N  =  yn . 

The  fluctuation  of  the  electron  density  with  respect  to  N  is  equal 
to  n  -  N  in  the  turbulent  fluid,  and  equal  to  -N  in  the  non-turbulent  fluid. 
On  the  average  (statistically),  the  density  fluctuation  is 

An  =  (n  -  N)6  +  (-N)  (1-6)  =  n6  -  N.  (9) 

Let  n  a  n  +  n  ,  where  n '  is  the  turbulent  fluctuation  of  n  (in  the  con¬ 
ventional  sense).  Then,  the  fluctuation  of  electron  density  is 

An  =  (6  ->')  n  +  n'  6  ,  (10) 

If  the  fluid  is  entirely  turbulent,  i.e.,  6  =  y  »  1,  An  becomes  n ' , 

When  the  turbulent  flow  shows  the  intermittency  phenomenon,  according 
to  Eq.  (11),  additional  fluctuation  of  electron  density  will  occur.  Hence, 
in  Eq.  (6)  the  fluctuation  in  electron  density  consists  of  two  parts:  the 
turbulent  fluctuation  and  that  introduced  by  the  intermittency.  The 
appearance  of  the  mean  electron  density  in  the  fluctuation  term  should 
be  noted. 

The  correlation  function  for  the  electron -density  fluctuations  is 

<(An)1  (An)2  >  =  (y^  -  <  fi2  >)  <  ^  *2  >  <  nj  n'2  > 

(11) 

(under  the  assumption  that  the  two  functions  6  and  n'  are  not  correlated) 


dill 
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where  the  subscripts  i  and  2  refer  to  points  at  r^  and  r^,  respectively, 
Let  the  correlation  coefficient  y  be  defined  by 


<  A,  62  > 

yl*2 


(12) 


Then,  the  correlation  function  can  be  written 


i  as 


<(&n),  (An)2  >  r  yjy 


(1  -  X)  n j  n,  +  V  <:  «2  > 


=  (1  *\)  f  y  <  N  '  N2  >  .  (13) 

The  problem  is  to  find  the  correlation  function  for  the  N  '  -  fluctuations 
and  the  correlation  coefficient  y. 

Fluctuations  in  dielectric  constant  may  arise  from  fluctuations  in 
temperature,  mass  density,  and  other  physical  properties  of  the  gas. 

But  for  the  wake  problem,  these  effects  are  believed  to  be  relatively 
unimportant  compared  to  that  given  in  Eq.  (11). 

(2)  Turbulent  fluctuations  ot  Electron  Density 

I  he  hypersonic  wake  flows  are  governed  by  the  processes  of 
mixing  between  the  inner  (turbulent)  wake  and  the  outer  inviseid  wake 
(see  reference  17,  tor  example)  so  that  large  gradients  of  the  mean 
velocity,  temperature,  mass  density  of  the  gas,  electron  density,  etc. 
in  the  direction  transverse  to  the  wake  axis  exist.  Considerable  progress 
has  been  made  recently  in  treating  this  problem  by  phenomenological 
turbulent-transport  theory  (references  17,  18,  19  and  20). 

The  basic  physical  process  for  the  production  of  electron-density 
fluctuations  is  assumed  as  follows:  Turbulent  motion  of  the  flow,  in  the 
presence  of  mean  electron-density  gradients,  will  generate  irregular 
fluctuations  of  electron  density.  Since  the  gradient  in  the  longitudinal 
direction  is  comparatively  very  small,  the  electron-density  fluctuation 
can  be  expressed  as  follows 
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N‘  ~  i  , 


(14) 


where  t  is  a  mixing  length  of  the  turbulent  velocity  field.  By  adopting 
von  (Carman's  idea,  this  length  is  assumed  to  depend  on  local  turbulent 
quantities:  a  turbulence  scale  (a  microscale  of  turbulence)  and  a 
turbulence  Reynolds  number  u^  X  tv  where  ur  is  the  component  of 
velocity  fluctuation  in  the  radial  direction  with  respect  to  the  wake  axis. 
In  particular,  if  the  relation 


l  -V 


(15) 


is  used,  the  expression  for  the  electron -density  fluctuation  becomes 

l2 

„<  “r  9N 

N  •  c  - - J7"  ■  04) 


where  c  is  an  absolute  constant  (of  order  1  presumbly  and  to  be  deter¬ 
mined  experimentally).  The  intensity  of  the  fluctuations  is 


where  C  is  another  absolute  constant,  and  q2  =  <  )u'2  +  v'2  4  w'2)  > 

is  the  intensity  of  the  velocity  fluctuations. 

In  considering  the  density  fluctuations  in  the  troposphere  and  the 

3 

ionosphere,  Villars  and  Weisskopf  took  i  as  a  scale  of  turbulence. 
Then,  the  intensity  of  the  electro-density  fluctuations  would  be  inde¬ 
pendent;  of  the  intensity  of  the  velocity  fluctuations.  This  result  may  be 
valid  for  the  ionosphere,  but  does  not;  appear  to  be  reasonable  for  the 
wake  problem.  Generally  when  turbulence  decays,  the  turbulence  scale 


10 


ncreases  but  it*  intensity  of  velocity  fluctuations  decreases,  from  the 

theory  of  decay  of  isotropic  turbulence  ,  in  the  initial  stage  of  decay, 

the  Reynolds  number  Xq/P  remains  constant.  Thus,  during  this  stage, 

i!  ^  X.  In  the  final  stage  of  decay,  q  -  t  "5^4,  and  X2  ~  t  (t  *  xUt 

1 

where  Uj  is  the  velocity  of  the  body).  Consequently,  in  the  final  stage, 

*■  3  /  4  ^  2  9 

Xq/P  -  t  .  The  presence  of  X  q  /V  in  Eq,  (17)  supports  the  idea 
that  the  intensity  of  electro -density  fluctuations  will  eventually  attenuate 
as  the  turbulent  energy  of  the  flow  decays.  In  Appendix  I,  another 
approach  is  used  in  attempting  to  justify  this  expression  for  the  intensity 
of  electron-density  fluctuations. 

In  the  above  consideration,  temperature  fluctuations  are  assumed 
to  have  no  direct  bearing  on  the  electron -density  fluctuations.  In  a  sense, 
the  turbulence  is  taken  as  frozen  chemically  as  far  as  the  electron-density 
fluctuations  are  concerned.  This  assumption  is  probably  reasonable  at  an 
altitude  around  150,000  ft,  in  standard  atmosphere,  at  which  the  recom¬ 
bination  rate  for  ions  and  electrons  is  very  slow. 

22 

Vi  liar  s  and  Weisskopf  also  looked  into  the  possibility  that  velocity 
fluctuations  will  produce  pressure  fluctuations,  which  in  turn  will  produce 
gas-density  fluctuations.  Assuming  that  the  electron-density  fluctuations 
follow  the  gas -density  fluctuations,  they  found 

?  2 

P  «  N  -S'  '  (18) 

VM 

where  u^  =  3p/p  .  The  intensity  of  fluctuations  given  by  Eq.  { 18)  will 
be  much  larger  than  that  given  by  Eq.  (17),  if 


1  hi s  condition  is  likely  valid  in  hypersonic  wakes,  since  v  X  tv  = 

M 

(qX/V)  (vM/q)  and,  for  hypersonic  wakes,  the  turbulent  Reynolds 

number  is  generally  much  larger  than  1,  and  v  /q  is  nearly  of  order  1 

M 
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Compared  to  an  expression  for  the  intensity  of  electron -density 
fluctuations  of  the  form  N*  m  N  (as  used  by  some  workers  for  under- 
dense  turbulent  scattering),  Eq.  (17)  shows  one  noteworthy  feature,  i,  e. 
it  contains  the  square  of  the  turbulent  Reynolds  number.  Therefore,  it 
is  a  sensitive  function  of  the  altitude  and  the  air  density.  It  is  proposed 
that  some  measurements  be  carried  out  to  check  this  point. 

(3)  The  Function  <  A.  A  > 

4  M 


In  Eq,  (13)  for  the  correlation  function  of  electron-density  fluctua¬ 


tions,  the  function  <(^6  >  appears.  This  is  the 
for  the  turbulence  probabilities  A  ,  and  6  at  r. 

4  w  J 


"c  orr elation  function  " 
and  r  ,  respectively, 


Very  little  is  known  at  this  time  about;  the  function  A  ,  but  experimentally 
<  Aj  A2> should  be  measurable  by  making  use  of  Townsend's  "A  "  signal 
method  *  . 


It  is  simpler  to  use  the  correlation  function  in  the  form  as  given 
by  Eq.  (13),  in  which  the  correlation  coefficient  \(R,  P)  appears.  At 
the  present  time,  however,  the  only  known  characteristics  about  x 
appear  to  be:  it  approaches  1  when  p  approaches  0  and  approaches  0 
when  0  is  very  large.  Tentatively,  the  following  form  may  be  used 


X  =  exp 


I4  +  T)  ‘ 


(20) 


where  I,  1}  and  (  are  the  coordinates  in  p  -  space,  R  (x)  is  the  radius 

2  w 

(or  half-width)  of  the  wake,  and  a  is  a  constant. 


(4)  Correlation  and  Spectrum  Functions 

In  a  previous  section,  the  following  expression  was  obtained  for 


the  electron-density  fluctuation 


N' 


>2 

A  ur  dN 
T - Vr" 


(16) 


The  question  is  whether  the  correlation  function  for  the  electron-density 

fluctuations  can  now  be  obtained  from  <  (u  )  (u  ),  >  .  Since  onlv 

r  1  r  2  ' 
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phenomenological  consideration  was  employed  in  deriving  Eq,  (16)  for  the 
excitation  mechanism,  of  electron-density  fluctuation,  this  relation  applies 
probably  only  to  the  "energy  containing  eddies"  of  the  N '  -  spectrum.  In 
the  high  wave-number  portion  of  the  spectrum,  dissipation  of  the  tur¬ 
bulent  energy  and  molecular  diffusivity  control  the  fine  structure  of  the 
turbulence  but  they  have  not  been  properly  taken  into  account  (see 
Appendix  I).  Therefore,  Eq,  (16)  will  be  used  only  to  find  the  intensity 
of  the  N  '  -  fluctuation,  as  given  in  Eq,  (17),  The  correlation  function 
for  the  N  -  fluctuations  will  be  assumed  to  be 


n.  n;  > 

i  L 


(R)  XN  (R,  0) 


(21) 


where  is  the  correlation  coefficient  for  the  N  -  fluctuations.  The 
spectrum  function  for  the  N*  -  fluctuations  miy  be  defined  by  the  Fourier 
integral 


i  k 


E  . .  (k ) 
N,  N 


Sir' 


<  N,  N,  >  e 

1  L 


dP 


(22) 


The  inverse  of  this  transform  is 


En  n  (it)  e  d  it 

.  m 

When  the  scale  of  the  electron-density  fluctuations  responsible  for 
the  scattering  (~  Xq/2  ,  where  X  is  the  incident;  electromagnetic  wave 
length)  is  in  the  high  wave -number  range,  results  from  turbulence  theory 
of  local  isotropy  may  be  used  for  the  correlation  and  spectrum  functions. 
It  is  necessary  to  take  y  =  1,  so  that  there  is  no  mtermittency  and  n '  a  N 
If,  in  addition,  there  exists  an  inertial  subrange  in  the  N'  -  spectrum,  the 
problem  becomes  more  tractable  and  has  been  treated  by  Batchelor23. 

The  relative  magnitude  of  the  molecular  diffusivity  D  and  the  kinematic 
viscosity  V  was  found  to  be  the  decisive  parameter  in  determining  the 


<N  N  > 

1  t 
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fofm  of  the  n  -  spectrum.  If  the  three -dim.en8io.otl  velocity  spectrum 
function  ia  of  the  form  as  follows  (Kolmogoroff  form) 

E  (k,  t)  -  <2/3  k'5/3  (23) 

then,  if  v  -  D,  in  the  same  k- range  the  corresponding  three-dimensional 

spectrum  function  E  ,  defined  by  E  *  4irk2E  ,  will  have  the  form 
n  n  n,n 

En(k<  ~  Cn  f  1/3  k'5/3  (24) 


In  the  above  expressions,  t  is  the  dissipation  of  the  turbulent  motion  by 
viscosity,  while  (  is  the  "dissipation"  of  the  n'  -  fluctuations  by  the 
molecular  diffusivity.  Both  these  two  functions  may  be  expressed  in 
terms  of  their  respective  intensity  and  scale  of  turbulent  fluctuations. 
The  scales  of  turbulence  for  the  velocity  and  n'  -  fluctuations  are  also 
of  the  same  order  (see  reference  10).  Under  this  circumstance,  the 
condition  required  for  the  presence  of  an  inertial  subrange  in  the  n'  - 
spectrum  is 


A± 

v 


>  >  1. 


(25) 


Since  the  size  A  of  the  "energy  containing"  eddies  can  be  assumed  to 
be  of  the  order  of  wake  radius,  by  using  the  relation  between  the  tur¬ 
bulent  Reynolds  numbers  (reference  10) 


Aq  A 
v  '  1* 


(26) 


an  estimate  for  the  value  of  Xq/v  in  the  turbulent  wake  can  be  made. 
The  precise  value  for  Xq/t'  at  which  an  inertial  subrange  will  be  in 
existence  is  not  known,  but  a  value  of  103  is  the  likely  minimum.  (The 
value  of  the  constant  A  in  the  above  equation  is  about  1). 

When  the  wave  length  X  ^  is  comparable  to  the  wake  width,  the 
larger  eddies  assume  the  role  of  scattering,  and  some  non-isotropic 
effects  of  turbulence  may  be  expected.  To  evaluate  these  effects,  the 


4 

» 


* 
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correlation  function  will  be  assumed  to  be.  of  a  form  which  alMr#  (the 
turbulence  scales  in  the  longitudinal  (along  the  wake  axis)  hhd  transverse 
directions  to  be  different.  This  idea 'is  adopted  from  Booker's  work  * 
on  radar  reflections  from  the  aurora.  Accordingly,  the  correlation 
function  for  the  N '  -  fluctuations  is  assumed  to  be  of  the  Gaussian  form 


< N|  Njj  >  (R)  exp 


4  ( 


jr 

L 


t72  + 

— 12 — 


) 


(27) 


where  I ,  fj  and  (  are  the  coordinates  in  p  -  space,  L  and  T  are  the 
correlation  lengths  in  the  longitudinal  and  transverse  directions, 
respectively,  and  N'2  ij  the  intensity  of  N '  -  fluctuations  at  R  » 

(1/2)  (r ^  +  r  ),  as  given  by  Eq.  (17).  The  Gaussian  correlation  function 
is  believed  to  be  valid  at  very  low  turbulence  Reynolds  numbers 
(corresponding  to  the  final  period  of  decay),  but  for  large  Reynolds 
numbers,  an  exponential  function  maybe  a  better  approximation  (see 
p.  173  of  reference  10).  There  is  some  uncertainty  as  to  the  values  of 

L  and  T  for  the  wake,  both  of  which  may  depend  on  R. 

25 

As  pointed  out  by  Tatarski,  use  of  spectrum  functions  for  the 
turbulent  scattering  problem  offers  Borne  advantage,  when  the  precise 
form  of  the  spectrum  function  is  known,  (in  the  inertial  subrange  for 
example).  The  important  question  of  how  the  scattering  cross  section 
depends  on  the  frequency  of  the  incident  electromagnetic  wave  can  be 


answered  without  any  ambiguity.  For  turbulence  at  smaller  turbulent 

Reynolds  numbers  (so  that  no  subinertial  range  exists), ''and  with  non- 

isotropic  effects,  the  exact  form  of  the  correlation  or  spectrum  function 

is  not  known.  Use  of  a  Gaussian  correlation  function  as  given  in 

Eq.  (27)  cannot  be  expected  to  predict,  in  general,  the  exact  frequency- 

dependence  of  the  scattering  cross  section. 

(5)  Some  Characteristics  of  the  Turbulent  Wake 

In  turbulent  wake  flows,  the  concept  of  self-preservation  as  intro- 
16  17 

duced  by  Townsend  '  has  been  found  to  be  particularly  useful.  For 
an  incompressible  axisymmetric  wake  flow,  the  mean  velocity  profile 


and  the  digtribution  of  the  turbulent  intensity  can  be  assumed  to  be,  of  the 
following  form 

U  *  Uj  -uQ  HtflQ)  ,  (28) 

2  "2  "  2  2 

u  4  V  4  w  .  uQ  g(r/f0)  ,  {29) 


where  Uj  is  the  free-stream  velocity,  r  *  (y2  4  x2)l/2,  and  u()  and 
*0*  suitably  chosen  scales  of  velocity  and  length,  are  functions  of  t 
only.  u0  is  called  the  defect  velocity,  when  the  function  f(r/*!0)  is 
chosen  to  be  unity  at  r  *  0,  (fQ  can  be  identified  as  the  wake  width). 

It **  found  that  self-preserving  solutions  are  admissible  by  the  equations 
of  motion,  and 


u 

0 

/  d  \2/3 

U!  ' 

V  54  -  *0/ 

l 

/x-x  \in 

*0 

(  0) 

1 

\  d  / 

where  d  is  a  typical  length  of  the  body  and  xQ  is  a  virtual  origin  of  the  flow. 

Wakes  generated  by  hypersonic  reentry  bodies  are  known  to  differ 
in  some  respects  from  their  low-speed  counterparts  .  Consideration 
must  be  given  to  a  turbulent  inner  wake  (when  transition  has  taken  place) 
and  an  inviscid  outer  wake.  Not  only  compressibility  effects  but  chemical 
reactions  involving  dissocation,  ionization,  recombination,  etc.  must  be 
taken  into  account  in  determining  the  mean  profiles  of  velocity,  tem¬ 
perature,  electron  density,  etc.  Progress  has  been  made  on  this  compli¬ 
cated  problem  (see  references  17,  19,  20,  for  example),  but  many  uncer¬ 
tainties  remain  to  be  resolved.  In  addition,  very  little  knowledge  is 
available  on  the  turbulence  structure  in  hypersonic  wake  flows.  Addi¬ 
tional  theoretical  and  experimental  studies  are  needed. 

As  far  as  the  scattering  phenomenon  is  concerned,  the  most  im¬ 
portant  physical  quantity  of  the  wake  is  the  electron -density  (mean) 


liilMiiiii 


'  "  . . .  . . . 


distribution.  According  to  Eq.  (17),  the  intensity  of  the  N*  -  fluctuations 
depends  on  the  gradient  of  this  distribution,  A  "wake  length"  dptihed  as 
the  distance  from  the  body  to  the  position  on  the  wake -axis  at  which  the 
plasma  frequency  is  equal  to  the  incident  electromagnetic -wave  frequency 
is  found  generally  to  be  a  useful  reference  length  (references  11  and  26), 
At  the  present  time,  uncertainties  in  the  electron-density  distribution 
and  the  wake  length  are  mainly  due  to  non -equilibrium  chemical  reactions, 
conditions  at  the  wake  neck  region,  turbulent  transport  properties,  and 
effect  of  Impurities  present  in  the  air.  It  is  Important  to  note,  however, 
that  by  referring  to  a  band  frequency  of  the  radar,  an  overdense  core 
surrounded  by  an  underdense  outer  portion  is  a  significant  feature  of  the 
wake.  A  scattering  model  with  this  feature  is  not  available  yet. 

There  are  strong  experimental  evidences  that  the  width  of  the 


turbulent  inner  wake  can  be  approximated  by  the  power  relation  (x  -  xQ)^^ 

as  in  the  incompressible  case  (see  references  17,  27  and  28).  If  the 

turbulence  scales  (microscales)  are  taken  as  constant  across  the  width 

of  the  wake,  and,  in  addition,  assumed  to  be  proportional  to  the  wake 

width,  then  both  L  and  T  in  Eq.  (27)  may  be  considered  as  proportional 
1  /  3 

to  (x  -  xq)  ,  This  would  be  in  variance  with  that  for  an  isotropic 

decaying  turbulence,  for  which  X~  t1^2  ~  (x-xQ)1^  .  However,  there 

is  some  experimental  evidence  that  for  axially  symmetric  wakes,  the 

i  n 

’’microscale  of  turbulences"  does  not  vary  according  to  x  (see 
reference  29). 

For  incompressible  wake  flows  (two-dimensional),  Townsend 
has  found  experimentally  that  the  turbulent  intensity  of  the  velocity 
fluctuations  in  the  turbulent  fluid  is  almost  uniform  across  the  wake, 
with  contributions  mainly  from  the  smaller  eddies.  This  result  should 
also  hold  true  for  hypersonic  wakes,  in  which  the  turbulence  mechanism 
is  expected  to  be  similar  to  the  low- speed  case  in  the  coordinate  system 
fixed  to  the  ground.  Hence,  according  to  Eq.  (17),  (in  which  X  is 
replaced  by  T  to  allow  for  the  non -isotropic  effect  in  turbulence  scales), 
the  variation  of  the  intensity  of  the  N '  -  fluctuations  follows  that  of 
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{  N/c)r)  in  both  the  longitudinal  and  the  transverse  direction*  iin.  the 
wake*  Thl*  result  point*  to  the  importance  of  having  an  accurate 
determining  of  the  mean  electron-density  distribution  in  hypersonic  wake, 
(6)  Summary 

The  correlation  function  for  the  electron-density  fluctuations  ha* 
been  found  to  be  of  the  form 


cfAnJj  (An)2  >  a  (1  -  +  n'2  x  Vn 


where  N  denotes  the  mean  electron  density,  is  the  intensity  of  the 
electron-density  fluctuations  as  given  by  Eq.  (17),  y  is  the  correlation 
coefficient  for  t  he  flow  inter  mittency  ,  and  x  N  t  hat  for  the  electron  - 
density  fluctuations.  Subscripts  1  and  2  refer  to  points  at  r,(  and  r  . 
for  which  the  double  correlation  coefficients  are  defined.  Suggested 
forms  for  the  coefficients  X  and  xN  are  Eq.  (20)  forx,  and  Eq.  (24) 

(for  the  spectrum  function)  or  Eq.  (27)  for  XN*  It  may  be  added  that 
the  expressions  for  the  electron  fluctuation  and  the  correlation  function 
apply  both  to  underdense  and  overdense  plasmas. 

IV*  scattering  cross  section 

The  correlation  function  for  the  electron-density  fluctuations  given 
by  Eq.  (32)  may  be  submitted  into  Eq.  (6)  to  obtain  an  expression  for  the 
scattering  cross  section  in  a  form  more  susceptible  to  physical  inter¬ 
pretation.  It  is  convenient  to  separate  the  cross  section  into  two  parts: 

e  a  b  I33) 


where 


l-X(R.P)  N(51)N(?2)e 


i  2k  •  P 


d  V 
1  V0  ’ 


1 2  ic  *  ft 

N71  W«R  X«.»)XN(R,0)e  dV  , 

V 
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And  2R  =  rj  +  r^,  P:  -r^.  If  the  effect  of  the  flow  intermittency  if 
omitted,  \  becomes  unity,  and  0^  vanishes ,  Then  (T^  reduces  to 


.  f  - 

i  '  n 


i  Zk  •  p 


(R)  dV 


In  considering  the  characteristics  of  the  scattering  cross  section, 

of  particular  interest  are  its  dependence  on  the  wave  number  (or  frequency) 

of  the  incident  electromagnetic  wave,  its  dependence  on  the  aspect  angle  ft 

(Figure  1),  and  how  the  various  physical  properties  of  the  wake  affect  the 

cross  section.  Evidently,  it  is  necessary  to  know  the  mean  electron- 

density  distribution  N  (r),  the  distribution  of  intensity  of  the  velocity 
2 

fluctuations  q  ,  the  turbulence  scales  L  and  T,  and  the  correlation 
coefficients  x  and  .  At  the  present  time,  none  of  these  quantities 
can  be  predicted  with  certainty.  Therefore,  the  results  obtained  in  the 
following  are  intended  to  show'  only  qualitatively  the  essential  features 
of  the  phenomenon  of  scattering  by  turbulent  wakes. 

(1)  Dependence  on  the  Wave  Number  and  the  Aspect  Angle 
By  referring  to  Figure  1,  write 


k  s  k(i  cos  ft  +  j  sin  ft)  (37) 

where  i  and  j  are  unit  vectors  in  x  and  y  (or  I  and  If)  directions, 

respectively.  If  the  volume  of  integration  in  Eqs.  (34)  and  (35)  is  very 

large,  the  dependence  of  the  cross  section  on  k  and  ft  can  be  found  by 

using  Fourier  integrals.  Since  the  functions  N,  \  and  attenuate  at 

^  50 

large  P,  from  the  theory  of  Fourier  integrals  (see  Titchmarch, '  for 
example)  the  scattering  cross  section  must  necessarily  decrease  as  the 
wave  number  k  (at  least,  for  sufficiently  large  values  of  R)  of  the 
incident;  electromagnetic  wave  is  increased.  In  addition,  if  the  attenua¬ 
tion  rate  is  not  isotropic  in  the  P-  space,  the  cross  section  will  become 
aspe c t; -an  gle  dependent . 


. . . . . 
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The  difference  between  0  ^  end  0^  in  their  dependence  on  k  and  0 


is  to  be  noted.  Conaider  firat  the  following  integral  in  ff 


a 


dV, 


i  2k  .  P 


Ntr^Nf^)  e 


dV  » 


12k(xcoi  0 +y*in0) 
N(x,y,x)e  dxdydss 


(38) 


Since  N  is  expected  to  attenuate  strongly  In  y,  but  weakly  in  x,  the  above 

integral  will  atemiate  strongly  in  k  sin  0  ,  and  weakly  in  k  cos  8.  Thus, 

cr  will  depend  very  strongly  on  k  for  8  ~(l/2)f,  but  will  be  leas  sensitive 

to  a  change  in  k  at  small  aspect  angles.  On  the  other  hand,  in  o  the  main 

contribution  to  its  k  and  0  dependence  is  due  to  the  difference  in  turbulence 

scales  between  the  transverse  and  axial  directions.  Since  this  difference 

in  scales  is  minor  compared  to  that  for  the  attenuation  rates  of  N  in  x  and 

y  directions,  the  aspect-angle  dependence  of  or^  will  be  less  prominent. 

The  simplest  case  is  obtained,  when  the  turbulence  is  isotropic  and  the 

three-dimensional  spectrum  function  is  of  the  form  given  by  Eq,  (24),  i. e. , 
-5/3 

E  k  ,  and  the  flow  intermittency  is  omitted.  Then,  from  Eqs,  (22) 

and  (36)  and  the  relation  E  *  4irk^E  , 

n  n,n 


ff  -  k 


-11/3 


(39) 


There  will  be  no  dependence  of  0  on  the  aspect  angle,  if  the  turbulence  is 
isotropic  the  scattering  volume  is  infinitely  large. 

Based  on  the  expressions  (33),  (34)  and  (35),  the  dependence  of  the 
scattering  cross  section  on  k  and  B  is,  in  general,  too  complicated  to  be 
expressed  in  a  simple  manner.  Further  analysis,  both  analytically  and 
numerically  is  needed. 

(2)  Some  General  Characteristics 

One  of  the  main  points  introduced  in  the  present  study  is  the  inter¬ 
mittency  phenomenon  of  turbulent  flows.  Since  the  correlation  coefficient 
Xfor  the  intermittency  is  always  smaller  than  1,  the  effect  of  the  flow 
intermittency  is: 
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(a)  to  introduce  a  contribution  or^,  which  depend®  directly 
on  the  mean  electron -density  distribution  N, 

(b)  to  reduce  the  contribution  from  the  N '  -  fluctuation*,  as 
seen  by  comparing  0fa  to  a"*1. 

Because  of  the  lack  of  knowledge  about  X  and  the  mathematical  difficulty 

in  evaluating  the  integrals  in  Eq.  (34),  it  is  not  possible  to  determine,  in 

general,  the  comparative  magnitude  between  a*  and  cr  +  (r,  ,  and  also 

h  b 

that  between  and  0^.  A  simple  analysis  by  using  mean  electron- 
density  distributions  of  an  exponential  type  is  given  in  Appendix  II, 

The  expression  (38),  multiplied  by  r  2  ,  would  appear  to  be  a 

3 1  ® 

coherent"  scattering  cross  section.  However,  the  contribution  or 

a 

arises  solely  from  the  flow  intermittency,  and  vanishes  if  \  is  taken  to 
be  unity.  In  addition,  N  is  the  mean  electron  density  averaged  over  both 
the  turbulent  and  non -turbulent  fluids,  i.e,  ,  the  intermittency  phenomenon 
is  also  included  in  N. 

For  the  wake  problem,  a  more  explicit  form  of  the  scattering  cross 
section  0  =  0  without  intermittency  (y  si,  n  »  N  and  n'*N') 

can  be  obtained  by  using  Eq.  (17)  for  n'“  and  Eq.  (27)  for  the  correlation 
coefficient  X  .  I  hen,  from  Eq,  (36),  if  the  volume  of  integration  is 
extended  to  infinity, 


0 


.1* 

_ 

i* 

n'  dVR 

exp 

V 

V 

cos  O  f  r?  sin  0) 


dfdijdt 


C^Zir) 


3/2 


2 

/  dn  \  t  _2 
{—)  L1  exP 

Sir**  2 

2 

2  A  ,  T2 

os  u  +  T  s 

in2  0) 

V 

X 

0 

dV„  , 


R 


(40) 


where  C  *  C  iV/j/2 


is  a  constant  (see  p.  10  ). 


Let 


n  (x,  :r)  =  F(x)  G(s) 


(41) 
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where  nQ  F(x)  i«  the  electron -density  diitribution  along  the  jivkfcf  axl#, 

•  ■  Cr/d2^3  (x  -  x  )  3],  both  F  and  G  are  non-dimensional,  and  G{o) !»  1. 
0 

w 

The  final  form  of  0  may  be  written  as  follows: 


®*  .  C  L  L  T  2.  2|  I  ,  <«> 

1  w  o  o  o  12 

where  is  the  wake  length  defined  in  a  proper  way,  and  are  the 
typical  scales  of  turbulence  in  the  axial  and  transverse  directions, 
respectively,  and 


Ij  a  2ir  G  (s)  s  d  s 


x  /L 
0  w 


r«,  .  LT2  8ff2  .  2  2  2  .  2  fl 

F(x)  - y  exp  -  - j  (L  cos  e  +  T  sin  o)  dx 

j  £» 


Ij  is  a  pure  constant,  while  1^  is  non-dimensional  but  depends  on  0,  L  /X  , 

T  A  ,  and  L  /x  . 
oo  wo 

^  „  _  2  4  2,2..  * 

Since  Cj  :  C  T  (T  q  IV  ),  the  cross  section  0  given  by 

4 

Eq.  (41)  is  proportional  to  ,  and  is  also  sensitive  to  the  turbulence 
Reynolds  number  T  q/v .  The  dependence  on  the  aspect  angle  0  vanishes 
when  the  turbulence  scales  L  and  T  are  the  same. 

To  be  able  to  evaluate  the  results  obtained  above,  it  is  necessary  to 
determine  the  quantities:  -  the  wave  length,  n  F(x)  -  the  electron 

density  distribution  along  the  wake  axis,  L  and  T  -  the  turbulence  scales, 
T  qfv  -  the  turbulence  Reynolds  number,  I  and  I  ,  and  the  proportional 
constant;  C.  How  these  quantities  vary  with  altitude  and  with  the  velocity 
and  geometry  of  the  body  will  be  of  particular  interest.  At  the  present 
time,  there  exists  considerable  uncertainty  in  the  electron -density  dis¬ 
tributions  in  the  wake.  Results  for  chemically  reacting  wakes  obtained 

19 

by  Bloom  and  Steiger  show  that  the  electron  density  is  not  sensitive  to 
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changes  In  altitude.  Lee*  and  Hromas  gave  the  results,  based  on 
frozen  and  equilibrium  conditions,  that  a  change  in  altitude  from  200, 000  ft. 
to  100,000  ft.  would  yield  electron  density  higher  by  about  two  orders  of 
magnitude,  Evidently,  much  work  in  this  area  is  needed,  before  mean¬ 
ingful  correlation  with  experimental  measurements  can  be  carried  out. 

(3)  Effect  of  Finite  Wake  Width 


Now  the  influence  of  finite  width  of  the  wake  on  the  aspect- 
dependence  will  be  treated  in  an  approximate  manner,  with  intermittency 
omitted.  Again  consider 


* 
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,*2  2  *2. 

ff 
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Assume  the  component  of  the  correlation  distance  in  the  transverse 
direction  is  limited  by  the  width  of  the  wake.  Then  the  integral  with 
respect  to  P  becomes 


1 


(2ff)  L  exp 


1 

A 

• 

w 

2  J. 

-8ir2(L2/A  2)  cos 2  0 

exp 

.  i  il 

0 

• 

21’2  K  o 

R 


dtjdC  . 


w 


where  R  is  the  radius  of  the  wake.  This  integral  can  be  evaluated 
w 


approximately  to  yield 
1 


4(2)  2  it  L  T2  exp 


o  2  ,,  2  a  la 

8ff  (L  /A  )  cos  B 
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erf  (R  /2  2  T)  *  I 
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where 


I  =  erf 


-  32  ff2  (sin2  9)(-J— )  {erf 
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2T' 


2  2  T 
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When  0  is  small,  higher  order  terms  in  sin  0,  etc.  can  be  neglected  in  I. 
The  above  result  indicates  that  the  back-scattering  cross  section  will  be 
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reduced  when  the  wake  width  is  finite.  In  addition,  lince  I  decreWf  f  ^hen 
0  is  increased,  the  effect  of  finite  wake  width  introduces  aspect -angle  ill 
dependence  into  the  scattering  cross  section.  The  fact  that  reduced  ttpis 
section  for  finite  wake  width  occurs  with  increasing  aspect  angle  seems 
to  be  reasonable  on  an  intuitive  argument,  namely,  the  Incident  electro¬ 
magnetic  wave  will  be  able  to  interact  with  more  eddies  at  small  incidence 
angles  than  at  larger  ones.  This  is  another  aspect  of  the  turbulent 
scattering.  In  the  latter  case,  the  scattering  cross  section  is  proportional 
linearly  to  the  number  of  electrons  anti  independent  of  the  aspect  angle  and 
the  wave  frequency.  In  the  above  consideration,  the  beam  width  of  the 
incident  wave  is  assumed  to  be  infinite. 

When  R  “*  »  ,  the  integral  becomes 


,3/2  2 

(2ff)  '  L  T  exp 

which  (the  same  as  given  in  Eq.  (40))  has  been  given  by  Booker  in 
reference  24,  This  result  written  in  the  form 

3/2  2 

(2ir)  '  L  T  exp 

shows  the  aspect -angle  dependence  of  non-isotropic  scattering  by  a  large 
volume.  If  L  >  T,  the  cross  section  increases  with  the  aspect  angle  in 
cont  rast:  to  the  case  of  a  wake  with  finite  width.  Therefore,  the  effect 
of  finite  width  is  found  to  be  quite  significant. 

V.  CONCLUSIONS  and  recommendation  for  further  work 

The  main  results  developed  in  this  study  are: 

(1)  The  phenomenon  of  flow  intermittency,  as  generally  recognized 
to  be  a  prominent  feature  of  turbulent  wake  flows,  will  introduce  electron- 
density  fluctuations,  in  addition  to  those  caused  by  turbulent  fluctuations. 


Sir  ,  2  2  2 
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(2)  The  electron-density  fluctuations  caused  by  turbulent  motions 

are  shown  to  depend  on  the  gradient  of  the  mean  electron-densky  distBibution, 
ft  turbulent  Reynolds  number,  and  a  turbulence  scale, 

(3)  By  using  the  above  two  results,  an  expression  for  the  bacfct 
scattering  cross  section  for  under  dense  wakes  is  given.  Two  contributions 
to  this  cross  section  arere2aa  and  re2ab  ,  given  by  Eqs.  (34)  and  (35)  , 

A  significant  point  is  that  depends  on  the  turbulent  Reynolds  number  but 
°a  does  not*  Thus,  at  higher  altitudes  aa  may  be  expected  to  be  dominant, 
while  at  lower  altitudes  ab  will  become  more  important. 

The  phenomenon  of  intermittency  is  found  to  introduce  significant 


aspect-angle  dependence  into  the  scattering  cross  section.  Other  effects 
discussed  in  some  detail  are  those  due  to  finite  wake  width  and  the  differ¬ 
ence  in  turbulence  scales  in  the  axial  and  transverse  directions  in  the  wake. 

Mathematical  difficulties  and  lack  of  knowledge  about  the  phenomenon 
of  the  flow  intermittency  make  a  thorough  analysis  of  these  effects  not 
feasible  at  this  time.  It  is  recommended  that 

(1)  further  theoretical  and  numerical  analyses  of  the  scattering 
problem  initiated  in  this  study  be  continued  in  order  to  develop  the  result 
obtained  so  far  into  more  useable  form,  and  to  correlate  them  with  experi¬ 
mental  work, 

(2)  to  study,  probably  by  experimental  means,  the  phenomenon  of 
flow  intermittency,  e.g. ,  to  measure  more  accurately  the  intermittency 
factor  for  hypersonic  wakes,  its  dependence  on  the  angle  of  attack  of  the 
hypersonic  body,  and  the  correlation  coefficient  for  the  flow  intermittency, 

(3)  to  undertake  basic  studies  of  the  turbulence  problem,  e.g. ,  the 
turbulence  theory,  turbulent  transport  properties,  intensity  of  turbulent 
velocity  fluctuations,  scales  of  turbulence,  and  other  turbulent  statistical 
properties  of  hypersonic  wakes,  and 

(4)  to  carry  out  an  analysis  of  turbulent  scattering  for  hypersonic 
wakes  based  on  the  model  that  an  overdense  core  is  surrounded  by  an 


underdense  outer  wake. 
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APPENDIX  I 


In  Section  III,  2,  the  following  expression  for  the  intensity  of  the 
turbulent  electron -density  fluctuations  is  given 


U7) 


Now,  a  partial  justification  is  attempted  here.  For  the  consideration  of 
the  excitation  of  the  electron -density  fluctuation  by  the  turbulent  velocity 
field,  effect  of  the  gradient  of  the  mean  flow  velocity  is  considered  as 
small  compared  to  that  of  the  mean  electron -density  distribution.  Assume 
the  electron  density  distribution  is  controlled  by  the  diffusion  process,  and 
use  a  coordinate  system  which  follows  the  local  mean  motion  of  the  fluid 
particle.  The  equation  governing  the  fluctuation  of  electron  density  is 


3N '  ,  dN 

t  ,uiie 


d2N' 


1  i 


3N'  „fl  3N' 
Ui  Tx”~  ‘  Ui  15x7" 


(A- 1) 


where  N'(x.,  t)  is  the  fluctuation  of  the  electron  density,  N(x.)  is  the 

mean  electron  density,  u.  is  the  velocity  fluctuation  and  D  is  the  molecular 

diffusion  coefficient.  The  second  term  at  the  right-hand  side  of  the  above 

equation  will  be  considered  as  small  and  omitted.  Similar  consideration 

has  been  developed  in  reference  32. 

Consider  an  initial  value  problem:  Let  N'  *  0,  at  t  *  0.  Coven  u 
(  i 

and  N,  determine  N'as  governed  by  the  differential  equation  (A-l). 

Solution  of  N  at  large  t  will  be  the  desired  electron-density  fluctuation, 

since  the  excitation  term  u.  ctN/3x,  and  the  dissipation  by  diffusion  will 

balance  each  other.  However,  since  little  is  known  about  u.  (x, ,  t)  and  the 

i  x 

analytical  behavior  of  N'  ,  the  equation  (A-l)  is  replaced  by 


(A-2) 


where  S  represents  the  excitation  term  and  X  is  a  turbulence  scale,  both 
assumed  to  be  time -independent  for  the  present  analysis.  The  solution 
of  Eq.  (A-2)  is 
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exp  (  -  X  t/D) 


(A- 3) 


Thus,  at  large  t. 


*  i  3N 

d  ir 


i  3N  V 


V  lx“  D 
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which  is  of  the  same  form  as  Eq.  (16).  The  constant  C  in  Eq,  (16}  is  thus 
probably  of  the  order  vfD,  which  is  independent  of  the  gas  density.  The 
intensity  of  the  fluctuations  will  be  essentially  the  same  as  given  by 
Eq.  (17)  for  the  wake  problem. 

If  Eq,  (A-l)  is  multiplied  by  N'  and  the  average  is  taken,  the 
following  equation  is  obtained  (only  significant  terms  are  kept) 
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(A-5) 


Now  introduce  the  eddy  diffusivity  Df],  by  the  relation 
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Eq.  (A-5)  becomes 
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Assume  Eq,  (A-7)  can  be  treated  in  a  manner  similar  to  Eq.  (A-l).  The 


result  is 
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Thus,  the  intensity  of  electron -density  fluctuations  Is,  approximately, 
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Comparison  of  Eq.  (A-9)  with  Eq.  (17)  shows  that  the  eddy  diffuslvity  is 
of  the  form 

2  ,2 

„  "  ^  X 

Dx  *  — jj - 5“  >  (A- 10) 

which  is  analogous  to  von  Harman's  concept  of  eddy  viscosity.21 

Following  Prandtl's  idea,  for  wake  flows  the  eddy  viscosity 

1 7 

coefficient  is  generally  assumed  to  be  given  by 


V  ~  (AU)  •  l 

i  o 


(A  - 11) 


where  AU  is  the  greatest  difference  in  the  mean  velocity  and  I  is  the 
width  of  the  wake.  Compared  to  the  thermodynamic  transport  coefficients 
D  and  V ,  the  difference  in  form  between  and  ,  as  given  by  (A- 10) 
and  (A -11),  respectively,  is  remarkable.  However,  the  only  justification 
usually  given  for  accepting  (A-ll)  is  that  satisfactory  mean  velocity  dis¬ 
tributions  in  the  mixing  zone  can  be  obtained  by  using  V T  of  this  form. 

It  remains  to  be  seen,  of  course,  whether  in  the  form  (A- 10)  will 
lead  to  satisfactory  results  for  the  electron-density  fluctuations.  A  point 
of  interest  is  that  DT  given  by  (A-10)  is  a  sensitive  function  of  altitude, 
since  V  depends  linearly  on  the  ambient  gas  density  while  f/D  is 
nearly  a  constant. 
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If  all  the  scales  L,  T,  a,  b,  P  and  Q  are  of  the  same  order  of  magnitude, 
ff  is  probably  larger  than  a  because  C  is  of  order  1  and  (  Tq/*)  >  1, 

*  A 

If  P  and  Q  are  very  small,  then  0.  /o  will  be  small. 

d  a 

The  more  significant  difference  between  0  and  a,  is  that  a, 

a  b  b 

depends  on  the  turbulent  Reynolds  number  T  q/v,  while  o  does  not. 

& 

Thus,  at  high  altitudes  <T  may  be  the  dominant  contribution,  and  at 

4 

relative  lower  altitudes  ctu  becomes  the  dominant  one. 
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The  or . HI  analytical  study  of  the  turbulent  scalier  tag  of  electro* 

ran  gnetic  wave*  is  directed  to  certain  feature*  of  the  phenomenon 
considered  as  peculiar  to  hypersonic  wake*.  Effect*  of  flow  intermittently, 
non-isotropy  of  the  turbulence  structure,  and  finite  width  of  the  wake  on 
the  scattering  cross  section  in  its  frequency  and  aspect-angle  dependence 
are  of  primary  concern. 

Phenomena  of  turbulence  revelant  to  the  scattering  problem  are 
briefly  diecussed.  In  particular,  it  is  shown  that  electron-density  fluctua¬ 
tions  of  a  turbulent  nature,  in  addition  to  those  caused  by  "turbulence", 
will  be  produced  by  the  interm ktlency  behavior  of  turbulent  wake  flows. 

This  fluctuation  Is  found  to  depend  on  the  mean  electron -density  distribution, 
and  Townsend's  intermittency  "ft"  function. 

Application  of  the  above  considerations  to  turbulent  scattering  by 
underdense  wakes  Ims  been  made.  In  contrast  to  the  Ionospheric  scattering, 
two  contributions  to  the  scattering  cross  section  are  obtained  in  the  present 
case:  the  first  one  arise*  from  the  intermittency  phenomenon  and  vanishes 
If  the  intermittency  is  omitted,  and  the  second  one  is  due  to  "turbulent" 
fluctuations.  The  frequency-dependence  of  these  two  contributions  are  seen 
to  be  quite  different.  It  is  further  shown  that  flow  intermittency,  non-isotropy 
of  turbulence  structure,  and  finite  width  of  the  wake  will  introduce  aspect- 
angle  dependence  into  the  scattering  cross  section.  A  possible  form  of  the 
total  scattering  cross  section  for  hypersonic  wakes  1*  also  given, 
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